The Paired Electron Crystal: order from frustration in the quarter-filled band 
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We present a study of the effects of simultaneous charge- and spin-frustration on the two- 
dimensional strongly correlated quarter-filled band on an anisotropic triangular lattice. The broken- 
symmetry states that dominate in the weakly frustrated region near the rectangular lattice limit are 
the well known antiferromagnetic state with in-phase lattice dimerization along one direction, and 
the Wigner crystal state with the checkerboard charge order. For moderate to strong frustration, 
however, the dominant phase is a novel spin-singlet paired-electron crystal (PEC), consisting of 
pairs of charge- rich sites separated by pairs of charge-poor sites. The PEC, with coexisting charge- 
order and spin-gap in two dimension, is the quarter-filled band equivalent of the valence bond solid 
(VBS) that can appear in the frustrated half-filled band within antiferromagnetic spin Hamiltonians. 
We discuss the phase diagram as a function of on-site and intersite Coulomb interactions as well as 
electron-phonon coupling strength. We speculate that the spin-bonded pairs of the PEC can become 
mobile for even stronger frustration, giving rise to a paired-electron liquid. We discuss the impli- 
cations of the PEC concept for understanding several classes of quarter-filled band materials that 
display unconventional superconductivity, focusing in particular on organic charge transfer solids. 
Our work points out the need to go beyond quantum spin liquid (QSL) concepts for highly frustrated 
organic charge-transfer solids such as k-(BEDT-TTF) 2 Cu 2 (CN) 3 and EtMe 3 Sb[Pd(dmit) 2 ]2, which 
we believe show frustration-induced charge disproportionation at low temperatures. We discuss 
possible application to layered cobaltates and |-filled band spinels. 

PACS numbers: 71.10.Fd,75.10.Kt,74.20.Mn 
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INTRODUCTION 



Strong Coulomb electron-electron (e-e) interactions 
can drive transitions from metallic to exotic insulating 
states, the most well known of which are the Mott- 
Hubbard semiconductor (MHS) and the Wigner crystal 
(WC). The MHS is a characteristic of systems with car- 
rier concentration per site p = 1 and is driven by strong 
onsite e-e repulsion, the Hubbard U interaction. Depend- 
ing upon the lattice structure the critical U at which the 
metal-insulator (MI) transition occurs, U c , can be + or 
finite^—. In contrast, the WC occurs in systems with 
p ^ 1, and is characterized by charge-ordering (CO) , 
i.e., a periodic arrangement of single charge carriers on 
the lattice. The WC is driven by strong onsite as well 
as inter-site Coulomb interactions^. Although in prin- 
ciple the WC is likely at any arbitrary p, it has been 
been studied most intensively for p = h bipartite lat- 
tices, where the nearest neighbor (n.n) Coulomb repul- 
sion can drive the MI transition^. The combined effects 
of e-e and electron-phonon (e-p) interactions are also of 
interest, usually in one dimension (ID), where the MHS 
can further exhibit the spin-Peierls (SP) transition. Im- 
portantly for our purpose here the above semiconducting 
states have been intensively studied over the past sev- 
eral decades, and are largely understood, although ar- 
guments regarding the magnitude of U c for formation of 
the MHS state in specific lattices or the detailed mech- 
anism of the MI transition may continue to persist. In 
the present work, we discuss a new correlated-electron 



semiconductor, the paired-electron crystal (PEC), that 
occurs in p = | systems in the presence of moderate 
to strong geometric lattice frustration 7 .. We believe that 
our work has direct application to 2:1 cationic or 1:2 an- 
ionic charge-transfer solids (CTS) that exhibit correlated 
insulator-SC transitions, and further applies to other in- 
organic strongly-correlated ^-filled materials. 

The combined effects of e-e interactions and geomet- 
ric lattice frustration are of strong current interest^. 
The bulk of the work here is for p = 1, where the Hub- 
bard model in the limit U — > oo reduces to the Heisen- 
berg spin Hamiltonian. Interest in the consequences of 
lattice frustration stems from the seminal proposal by 
Anderson that the ground state of the Heisenberg an- 
tiferromagnet (HAF) model on a triangular lattice is a 
quantum spin liquid (QSL) with no spin ordering even 
at zero temperatur o 10 ' 11 . The type of wavefunction usu- 
ally assumed to describe a QSL is often referred to 
as a resonating valence bond (RVB) state. Whether 
or not RVB states appear in the square lattice for p 
slightly different from 1, and the relationship of such 
states to superconductivity (SC) in doped strongly corre- 
lated semiconductors remains contentious. An extension 
of the RVB theory of dopant-induced SC in p ^ 1 is 
the proposal that frustration-induced SC occurs in the 
anisotropic triangular lattice within the simple Hubbard 
model even for p exactly 1 , where a narrow superconduct- 
ing phase is straddled on both sides by broader param- 
agnetic metallic (PM) and antiferromagnetic (AFM) in- 
sulator phases^—. It has been claimed that this transi- 
tion explains the SC in the CTSi2~— . Recent numerical 
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work by us and others, however, have determined that 
SC is absent within the p = 1 triangular lattice Hub- 
bard model^ir— and the earlier results are artifacts of 
mean-field approximations. 

While the ground state of the HAF on the isotropic tri- 
angular lattice is now known to be the ordered 120° AFM 
rather than the originally proposed QSL state^Zr— , other 
frustrated lattices, most notably the Kagome lattice^, 
have been investigated in the search for QSL states. 
Proposed ground states here include various types of 
QSL states^— as well as valence-bond solid (VBS) 
states^— . The literature on VBS states has a long 
history going back to the well known Ghosh-Majumdar 
model^. The common theme in works on VBS is the 
frustration-driven transition from the AFM state to a 
total spin S = singlet state. We have found a similar 
frustration-driven transition from the AFM to a S — 
state in strongly correlated systems with p = i, where 
reduction to a spin Hamiltonian is not possible. 

In contrast to the voluminous literature on correlated 
and frustrated systems at p = 1, the literature on frus- 
trated p 7^ 1 is relatively sparse and new. The discovery 
of SC— in hydrated Na^CoO^ has spurred interest in 
correlated systems away from p = li2r— , although to 
the best of our knowledge only isotropic triangular lat- 
tices have been studied. We will specifically focus on 
P = \ within the present work — on the triangular lattice 
with varying anisotropy. In the square lattice limit at 
this density, we show that spontaneous in-phase dimer- 
ization occurs in the presence of electron-phonon inter- 
actions modulating n.n hopping integrals, leading to an 
effective p = 1 system with one electron per dimer and 
AFM order. This result is the origin of the so-called the 
dimer Mott-Hubbard model 4 -^ 4 - that is commonly used 
to describe the 2:1 cationic or 1:2 anionic organic CTS. 
Very recently we have proposed that under the influ- 
ence of lattice frustration this dimer Mott-Hubbard AFM 
state gives way to a spin-paired state that we termed the 
PEC ' . The PEC is different from any of the above more 
well known correlated semiconducting states in that it 
is a WC of Heitler-London spin-singlets — simultaneously 
charged-ordered and spin S = 0. Alternately, the PEC is 
the p — h equivalent of the p = 1 VBS. A conceptually 
similar state was postulated for the electron gas many 
years back by Moulopoulos and Ashcroft 45,46 . There is 
a fundamental similarity between this earlier work and 
ours, in that in both cases the pairing is driven by the 
exchange interaction. Our previous work£ only consid- 
ered a limited set of parameters, and the full phase dia- 
gram was not discussed. Here we give a more complete 
phase diagram, including the competition with the WC 
state that was ignored before. As with the VBS state at 
p = 1, the PEC is a consequence of frustration- induced 
quantum effects. Its extraordinary stability at p = h 
is a commensurability effect (recall that MHS and WC 
formation also require commensurability). 

In section |TT] we introduce the model we consider. We 
include e-p interactions to stabilize the lattice dimeriza- 



tion that gives the effective p — 1 dimer lattice. We 
present the physical mechanism behind the PEC forma- 
tion by briefly discussing simple molecular clusters, for 
which we show that spin-singlet formation in p = ^ nec- 
essarily requires charge disproportionation. The charge 
disproportionation in the infinite one-dimensional (ID) 
chain and the so-called zigzag ladder leads to periodic 
CO, viz., the simplest PECs. Following the discussions 
of these simple cases, we introduce the two-dimensional 
(2D) lattice that will be the focus of this work, and 
discuss the different possible phases. In section [Hi] we 
present out numerical results for the 2D system, cov- 
ering a wide region of parameter space. In section IIVI 
we discuss the relevance of our results for several classes 
of p = g materials and the outlook for understand- 
ing unconventional SC. We particularly emphasize the 
cases of the organic charge-transfer solids CTS k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 and EtMe 3 Sb[Pd(dmit) 2 ]2 which have 
been described as QSLs within the effective p = 1 
scenario^ - — . We believe frustration-induced charge dis- 
proportionation is an alternate possibility. 

II. THEORETICAL MODEL 

A. Hamiltonian 

The Hamiltonian we consider contains electron hop- 
ping, semi-classical inter-site and onsite e-p couplings, 
and onsite and n.n. Coulomb interactions: 

H = - t v {l + a v A ij )B ij +^ Y. K^ijiX) 

i i 

• <y> 

In Eq. [TJ v indexes the different bond directions in the 
lattice; for example v = x in ID and v — {x, y} in the 
2D square lattice. Our actual calculations (see below) are 
for the anisotropic triangular lattice, v — {x,y,x + y}. 
Bij — ^2 a {c\ a Cj a + H.c.) is the electron hopping between 
sites i and j with electron creation (annihilation) opera- 
tors c\ a {ci a ). a v is the inter-site e-p coupling constant, 
K£ is the corresponding spring constant, and Ay is the 
distortion of the bond between sites i and j. Vi is the 
intra-site phonon coordinate and (3 is the intra-site e-p 
coupling with corresponding spring constant Kp. Both 
Aij and Vi are determined self-consistently— a v are 
in general taken close to the minimum value needed for 
the transition to occur, our goal being the replication of 
the same instability from finite cluster calculations that 
would occur in the infinite system for + coupling. U and 
Vtj are on-site and n.n. Coulomb interactions, respec- 
tively. The physically relevant range of is Vij < ¥ 
based on comparison between p = 1 and p = \ CTS^i. 
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FIG. 1: Charge difference An between members of the same 
dimer as a function of the hopping integral t' corresponding 
to the dotted bonds, for the two molecules (a) and (b) given 
as inserts. Each 4-atom molecule contains two electrons. The 
intra-dimer double bonds have strength t\ = 1.5, and the 
single bonds in (b) are £2 = 0.5. The results shown are for 
U = 4 and V = 0. Solid (dashed) curves show An for the 
linear (square) molecules. Lines are guides to the eye. Filled 
and empty circles of the molecules correspond to sites with 
charge densities 0.5 + An/2 and 0.5 — An/2, respectively. 

B. Coupled spin-singlet and CO at p = | 

We first present a simple qualitative discussion of cou- 
pled spin-singlet and CO formation at p = |. The ideas 
are quite general and we argue that the mechanism is 
independent of dimensionality. The key requirement is 
that the density must be exactly p = |, as the effect re- 
quires commensurability. Consider a single dimer of two 
sites with one electron. The electron populations per site 
are 0.5 each, but the quantum mechanical wavefunction 
for the system is the superposition -i=|10 + 01), where 
1 and are site charge densities. If one now brings two 
of these dimers together, as in insert (a) of Fig. [I] the 
composite wavefunction of the two-dimer system can be 
written as ||1010 + 1001 + 0110 + 0101). If the two elec- 
trons are in a spin-singlet state then within the simple 
Hubbard Hamiltonian the configuration 0110, in which 
singlet stabilization can occur from a single n.n. hop that 
creates a virtual double occupancy, must dominate over 
the configurations 1010 and 1001, in which singlet stabi- 
lization requires two and three hops, respectively. Thus 
as the singlet bond between the dimers gets stronger we 
expect a charge difference An between sites belonging to 
the same dimer (nominally between sites 1 and 2, or be- 
tween sites 3 and 4 in the linear chain of Fig. [IJ . While 
some charge disproportionation must occur in finite lin- 
ear chains from end effects alone, we note that our pro- 
posed picture demands that similar charge disproportion- 
ation occurs between members of the same dimer even in 
the case of the periodic molecule shown in the insert (b) 
of Fig. [1] In this case the charges on the sites connected 
by the diagonal bond must be larger than 0.5, while the 
charges on the two other sites must be smaller. Impor- 
tantly, the modulation of charge density, bond orders, as 
well as spin-singlet pairing all occur cooperatively, and 
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FIG. 2: (a) 2D lattice used for calculations in this paper, 
a square lattice with hopping t x = t y = t and frustrating 
bond t x + y = t' (dashed lines), (b) Dimerized lattice. Double 
(single) lines indicate stronger (weaker) bonds. Site charge 
densities are uniform in both (a) and (b). (c) The PEC state 
as it occurs in this lattice. Filled (open) circles correspond to 
sites with charge density p = 0.5 + 5 (p = 0.5 — 8)^. Heavy line 
shows the location of singlet-paired sites, (d) Wigner crystal 
charge ordering occurring for large V x and V y . (e) Wigner 
crystal-spin gap phase with bond alternation along diagonal 
directions. See section [ill Bl 



any of these observables may be used as an order param- 
eter in the case of a real transition. In Fig. Q] we have 
plotted An versus the hopping integral t' corresponding 
to the dotted bonds in the molecules shown in the in- 
sert, for the ground spin-singlet state. In both cases, as 
t' increased from zero An becomes nonzero and increases 
with t' . Importantly, in the spin-triplet S = 1 state the 
sign of An is reversed in the cyclic molecule, indicating 
repulsive interaction among the electrons. 

Coupled CO and singlet formation occur also in the 
SP state in ID p = i systems^—. A key difference 
from the SP transition in p = 1 is that for p = i, a 
MI transition first occurs at a intermediate temperature, 
followed by the SP transition at low temperature. For 
p = |, the MI leads to either a bond-dimerized or to 
a charge-ordered WC state with equal bond lengths 57 . 
Provided the n.n. Coulomb interaction is not too strong, 
a SP transition occurs from either insulating stated, re- 
sulting in a ground state with period-4 CO • • • 0110- • • . 
This state is the simplest realization of the PEC and 
may visualized as a second dimerization of dimer units 
of molecules; the singlet bond giving the spin gap (SG) 
forms between adjacent dimer units. The intermediate 
temperature bond-dimerized states, WC state, and PEC 

state with SG are all found experimentally in quasi-lD 
CTS5£ 

Zigzag ladder systems, coupled two-stack systems in 
which each site on one stack is coupled to two sites on 
the other stack, are a second realization of the PEC state. 
CTS zigzag ladder materials that are p = \ have been 
found with spin-gap transition temperatures much larger 
than in ID p = \ SP materials 5 ^. The insulating ground 
state in this case may be understood again as a PEC state 
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occurring in a zigzag ladder lattice, with singlet bonds 
oriented between the two chains^. Unlike the ID case, 
bond orders are now modulated in several lattice direc- 
tions, leading to a larger SG than for the ID PEC case 60 . 
Interestingly, in both the linear chain and the zigzag lad- 
der, the p = i PEC is obtained by simply removing al- 
ternate spin-singlet bonds from the corresponding p = 1 
VBS, and replacing them with pairs of vacancies. 



C. Competing broken symmetries in 2D 

In this paper we will focus on ground state solutions 
of Eq. [JJ in 2D in the presence of variable lattice frustra- 
tion. The lattice we choose is a 2D square lattice with 
a single frustrating bond, as shown in Fig. Ufa). Thus 
v = {x, y,x + y} within Eq. [TJ for this lattice. In most 
of the results we will present, t x = t y = t, although we 
will also consider t x ^ t y in some cases. Energies will 
be given in units of t. We will take the frustrating bond 
t x+y = t' in the range < t' < 1, covering the wide 
region between the unfrustrated square lattice (t' = 0) 
and the nearly isotropic triangular lattice (t' — 1). For 
the inter-site e-p coupling, unless denoted otherwise we 
choose a x = a y = a and a' = 0, with similarly identical 
spring constants — = K a . For all calculations we 
assume periodic boundary conditions. 

Our calculations are largely for V x — V y , but vari- 
able V x+y = V. In Reference we presented limited 
numerical results for a select set of Coulomb interaction 
parameters (U = 4, V x = V y = 1, V' = 0) demonstrat- 
ing transition from Neel antiferromagnetism to the PEC 
state in this lattice when t' exceeds a critical value t' c . 
For completeness and for giving an introduction to the 
various competing states in 2D we briefly review these 
results here. For small t', the self-consistent solution of 
Eq. [1] gives spontaneous dimerization along the x axis as 
shown in Fig. HJh) . The dimerized lattice state is effec- 
tively ^-filled with one carrier per dimer and has Neel 
AFM order between dimers for finite U . The Neel order 
is very clearly observable from the spin-spin correlations 
calculated for 4 x 4 clusters 7 -. Importantly, the charge 
density (rii) of all sites is exactly 0.5. 

As t' increases, frustration reduces the strength of 
the antiferromagnetic correlations. Provided the n.n. 
Coulomb interaction Vij is not too strong (see below), 
at a critical t' — t' c antiferromagnetic correlations dis- 
appear and charge disproportionation develops, with the 
charge densities within each dimer becoming inequiva- 
lent. This is shown in Fig.^c). In this state, the charge 
densities follow the pattern ■ • • 1100- ■ • along the x and 
x + y directions, and the pattern • • • 1010- • • along the y 
direction. The bond distortion is period-4 along x and 
period-2 along y. The strongest bond orders (By) oc- 
cur between adjacent charge-rich '1-1' sites in the x + y 
direction. Importantly, the spin-spin correlations change 
dramatically for t' > t' c . They are strongly negative be- 
tween the bonded '1-1' sites along the x + y direction (see 



Fig-HI c )), and are nearly zero between either member of 
the pair and all other sites, indicating the formation of 
spin-singlet bonds 7 . Any of the observables An, bond 
order between the charge-rich sites, or z-z spin correla- 
tions, (SfSj), may be used as order parameters for the 
PEC stated. 

In section IIIII we present further details of the PEC 
phase and the full parameter dependence of Eq. [TJ with 
the goal of demonstrating that (i) the transition to the 
PEC that we are interested in is driven by quantum ef- 
fects due to frustration only; and (ii) the PEC occurs 
over a broad region of parameter space. Given the num- 
ber of parameters in Eq. [TJ it should be relatively easy 
to generate CO driven by specific (presumably artificial) 
choices of Vij. Such classical results would be uninter- 
esting. We therefore consider several distinct choices of 
Coulomb interactions: (i) U > and all Vij = 0, (ii) 
U > 0, V x = Vy = V, and V = 0, (iii) U > and 
V x = V y = V = V. PEC formation occurs in all of these 
parameter regions. We also show that for sufficiently 
strong n.n. Coulomb interactions corresponding to pa- 
rameter region (ii), the WC phase with checkerboard CO 
(Fig. E^d)-(e)) is the ground state of Eq. [TJ We also will 
consider other modifications of the basic lattice, viz., sign 
oft' opposite site to t, and t x ^ t y . Finally, we will argue 
that our results are not consequences of finite size effects 
and are to be expected in the thermodynamic limit. 

III. RESULTS 

A. U > 0, V s = V v = V = 

Figs. EJa)-(d) show the charge disproportionation An, 
bond order (Bij), n.n. z-z spin-spin correlation {SfSj}, 
and spin gap A$t as a function of t' . A$t is defined 
as the excitation energy from the ground 5 = state to 
the lowest S = 1 state. U here is finite but all V terms 
are zero. Sites i and j in Figs. G3b)-(c) correspond to 
two T' sites in the PEC state connected by a t' bond 
(filled circles connected by heavy lines in Fig. [5Jc)). As 
with the results in Reference which included the V x and 
V y interactions, for small t! the charge density is 0.5 on 
all sites, and antiferromagnetic order can be seen in the 
spin-spin correlations (not shown here, see Fig. 3(a) in 
Reference 0) . Nonzero t' c in Fig. [3] is a consequence of 
the nature of the diagonal bonds in Fig. [UJb); the diag- 
onal bonds inside each plaquette with two strong dimer 
bonds actually strengthen the AFM, and only the inter- 
plaquette diagonal bonds have a frustrating effect. At 
small t! these two effects appear to cancel, and there is 
only a weak effect on the AFM. 

For t' > t' c , An becomes nonzero, with the charge 
pattern as shown in Fig. dfc). Similarly, the bond or- 
der between paired '1-1' sites increases abruptly, and 
the z-z spin-spin correlation between these sites becomes 
strongly negative, and nearly zero with all other lattice 
sites (see Reference 0) , indicating formation of a singlet 
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FIG. 3: (Color online) Order parameters for the 4x4 lattice versus t! . Parameters are a = 1.1, (3 = 0.1, and K a =Kp=2. Circles, 
squares, diamonds, and triangles are for U = 2, 3, 4, and 6, respectively. V x = V y = V' = for (a)-(d), and V x = V y — 1 and 
I 7 ' = for (e)-(h). (a) and (e) show the charge disproportionation An, (b) and (f) the bond order Bi t j between charge-rich sites 
i and j connected by t' (see Fig.f2Jc)), (c) and (g) spin-spin correlations between these sites, and (d) and (h) the singlet-triplet 
gap Ast- In all cases the CO pattern is as shown in Fig. [2jc). For all plots lines are guides to the eye. 



bond. Although Ast is nonzero in all cases in a finite 
cluster, we nevertheless see a large jump in Ast at t' c , 
also indicating spin-singlet formation. The increase in 
bond order strength and strength of spin-spin correlation 
clearly follow the same pattern as An. 

These results show that the n.n. Coulomb interaction 
is not essential for formation of the PEC state. Unlike 
the WC phase where CO is driven by the n.n. Coulomb 
interaction, the PEC state is a consequence of geometric 
lattice frustration. Increasing U moves t' c to larger t': 
U tends to strengthen the AFM phase and therefore in- 
creasing U is expected to make the AFM order persist for 
stronger lattice frustration. The transition also becomes 
more discontinuous to changes in t' as U increases, sug- 
gesting it may be continuous for small U and first order 
for large U. 

In both the ID and in the zigzag ladder lattice, the 
PEC state occurs unconditionally even in the noninter- 
acting limit (U = Vij = 0) for any finite e-p coupling. In 
these two cases, the unconditional occurrence of PECs 
is a consequence of simple nesting. In contrast, in the 
isotropic p — | 2D band considered here, the lack of 
nesting forbids an unconditional Peierls transition. In 



agreement with this, we found that for U < 1, the PEC 
phase did not occur. Instead, the self-consistent calcula- 
tions converged to disordered states with no clear charge 
pattern. This is an indication that in the thermodynamic 
limit, the preferred ground state is one of uniform charge. 
This result is reminiscent of that in quantum spin sys- 
tems: the simplest VBS transition, the SP transition, 
can be predicted from nesting behavior in the ID XY 
model following Jordan- Wigner transformation 61 . This 
is, however, not true in 2D frustrated spin systems. 



B. U > 0, V x = V y > 0, V' = 

We next consider the effect of n.n. Coulomb interac- 
tions V x = V y = V, but V = 0. Figs. GHe)-(h) show the 
same order parameters as in Figs. [3Ia)-(d) for V = 1. 
Comparing the data with and without V, the effect of 
moderate V is to strengthen the PEC state — the mag- 
nitude of all order parameters increase when V > for 
a fixed value of U. For fixed U, the AFM-PEC bound- 
ary t' c also moves to smaller t' with increasing V. Fig. 0] 
shows the phase diagram in the t'-U plane for both the 
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FIG. 4: Phase diagram as a function of t' and U, for j3 = 
0.1, and K a = = 2. Squares are the AFM-PEC phase 
boundary for a — 1.1 and V x = Vy = V = 0; Diamonds are 
for a — 1.1, V x = V y — 1, and V' = 0; Circles are for a — 1.2, 
V x = V y = 1, and V" = 0. For each case, filled (open) points 
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V = and V > cases. Here, for each value of U, if 
was varied until the AFM-PEC transition occurred, and 
the first t' where An became nonzero was taken as the 
AFM-PEC boundary. 

One expects that when V is above a critical value V c , 
the • ■ • 1100- • • PEC CO will give way to the checker- 
board WC state. In ID, this transition occurs exactly 
at V c = 2 in the limit U — > oo, and at a larger V c 
for finite U— Previous exact diagonalization for a 2D 
cluster- 2 (without however e-p interactions or t' as con- 
sidered here) found V c ~ 2.1 for U — 10, and showed that 
V c increases when V > 0. 

Fig. [5] shows the evolution of An and diagonal bond 
orders with V for U = 6. While difficult to see in Fig. [S] 
due to the choice of axis scales and parameters, in the 
PEC region An increases with increasing V. For the 
parameters of Fig. G2 (17=6, t' = 0.8, a = 1.1, = 0.1, 
K a = Kp = 2) the charge order pattern changes from 
the PEC (Fig. |2(c)) to the checkerboard WC (Fig. [2(d)- 
(e)) at V » 1.52. In addition, An increases sharply when 
entering the WC phase. At larger V ~ 2.0 a slight cusp 
occurs in the An versus V plot. At the same time, the 
pattern of bond orders changes: for V < 2 in the WC 
phase the bond orders alternate strong-weak along the 
x+y direction (as shown in Fig. [2(e)), while for V > 2 the 
bond orders along x + y are uniform (Fig. [2(d)). Within 
the WC phase region there are therefore two sub-phases: 
a phase with equal length bonds in the diagonal x + y 
(t') directions, and a phase in which these bonds become 
dimerized. The added bond dimerization in the diagonal 
direction will result in a spin gap, and we denote this 
phase as the Wigner Crystal-Spin gap (WC-SG) phase. 
A similar spin-gapped WC phase can be found in a small 
region of parameter space in the ID mode l 56 i 63 . 

Fig. [5(a) shows the resulting phase diagram in the t'- 

V plane for J7=6. The V c we find in the t' = limit is 
slightly larger (U = 6, V c s» 2.6) than the results of Ref- 
erence 62 (U = 10, V c w 2.1); however, both the smaller 
U here as well as the e-p coupling in Eq. [TJ would be ex- 
pected to increase V c . Unlike the AFM, PEC, and WC 



FIG. 5: Variation of order parameters with V = V x = V v 
{V = 0) for U = 6, t' = 0.8, a = 1.1, K a = Kp = 2, p = 
0.1. (a) charge disproportionation An (b) bond orders along 
two successive t' bonds. In PEC region, these correspond to 
bonds between sites with '1—1' (squares) and '1-0' (circles) 
occupancy. In the WC-SG and WC regions, both bonds are 
between '1—1' sites. Lines are guides to the eye. 
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FIG. 6: (a) Phase diagram for the 4x4 lattice as a func- 
tion of t' and V = V x = V y , V = 0, for U = 6, a = 1.1, 
P — 0.1, and K a = Kp = 2.0. (b) Same as (a), but with 
V — Vx ~ V y = V . For both (a) and (b), points between an- 
tifcrromagnetic and PEC phases are determined as discussed 
in section III Al The boundary between PEC, WC-SG, and 
WC phases is discussed in the text. Lines are guides to the 
eye. 



phases, the WC-SG phase is limited to a relatively nar- 
row range of parameters. Larger lattice calculations are 
needed to confirm whether the WC-SG phases persists in 
the thermodynamic limit. 



7 



C. U > 0, 14 = V y = V > 

The V interaction destabilizes the checkerboard- 
pattern WC, leading to a metallic phase in the absence 
of e-p interactions^. Here we will consider parameters 
V x = V y = V = V. For EqQ] without e-p interactions, 
U = 10, and t' < 0.1, exact diagonalization found— that 
in this case a metallic phase exists for V up to at least 
V —5. Charge fluctuations within the metallic phase 
adjacent to the WC were speculated to cause a CO-to- 
SC transition^. In our calculations, summarized in the 
phase diagram in Fig.[6jb) for U = 6, we also found that 
the WC phase does not occur, but rather than being 
metallic the system is insulating — either AFM at small 
frustration or PEC at large frustration. In this case the 
AFM-PEC transition can occur over a wide range of lat- 
tice frustration, 0.2 < t' c < 0.7. Within the PEC phase 
the CO pattern remains the same for all V, although An 
increases with V as in the V = case considered in the 
previous section. Although our calculations are for one 
value of V only, Fig. Hfb) suggests that the PEC region 
is broadened relative to that in Fig. [B]Ja) for any V ^ 0. 

D. Bandstructure and electron-phonon coupling 

As shown in the previous section, variation of Coulomb 
interactions can cause a substantial variation in the ex- 
tent of frustration needed to form the PEC state. Next 
we show the effect of varying the one-electron parameters 
in Eq. [TJ t v and a. v . 

Due to the lack of particle-hole symmetry in the 
anisotropic triangular lattice, differences might be ex- 
pected when t' is taken as negative. However, as Fig. [4] 
shows, we found only a small variation in the AFM-PEC 
phase diagram when the sign of t' is changed. This is con- 
sistent with the expected mechanism for spin frustration 
in an effectively ^-filled band: the frustrating exchange 
interaction is proportional to (t') 2 , so reversing the sign 
of t' should only change the effective frustration at higher 
order. 

Fig. S] also shows the effect of changing a. As ex- 
pected for a cooperative transition, stronger e-p cou- 
pling increases the size of the PEC region. The effect 
of anisotropy of the inter-site e-p interaction is shown in 
Fig. E^a), where a x is fixed at 1.1 and a y is varied. In- 
creasing either a x or a y separately strengthens the PEC, 
as shown in Fig. EJa), where the a y /a x ratio is varied. 
The AFM-PEC phase diagram is relatively insensitive to 
the ratio of t y /t x . In Fig. EKb), t x is fixed at 1 and the 
value of t y is varied — the resulting t' c does not depend 
sensitively on the choice of t y /t x . 

Summarizing sections iHlT aVf d) . we see the PEC state 
in our calculations for a wide range of parameters. In 
every case, we first arrived at parameters that placed the 
system in the AFM or WC phases in the t' = limit, 
and then varied only t' . Thus the transition to PEC 
is a consequence of frustration alone. With variation of 
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FIG. 7: Phase diagram variation on hopping and electron- 
phonon interaction anisotropy. In both (a) and (b), [7 = 6, 
V x = V y = 1, V = 0, = 0.1, and K a = Kp = 2. In (a), 
t x = t y — 1 and a x = 1.1, but a y is varied. In (b), t x — 1 and 
t y is varied, with identical a x — a y — 1.1. Lines are guides 
to the eye. 



e-p coupling and Coulomb interactions, the amount of 
frustration needed to drive the AFM-PEC transition can 
vary over a sizable range of frustration. We will discuss 
the implications of this further in relationship to the CTS 
materials further in section ITVl 



E. Finite-size issues 

We have performed several checks on our calculations 
that indicate that the PEC state found in our numeri- 
cal calculations is an intrinsic property of Eq. [1] and not 
induced by finite size effects. 

(i) Noninteracting bandstructure: One common finite- 
size effect in numerical calculations are changes in the 
Fermi-level degeneracy or level crossings of the nonin- 
teracting system. In the lattice of Fig. [2fa) the Fermi 
level degeneracy does not change throughout the range 
< t! < 1, remaining 2-fold degenerate throughout this 
range. This degeneracy is broken by ir-axis dimcrization, 
shown in Fig. EJb) , giving a nondegenerate Fermi level 
for < t 1 < 1. At t' = 1 the degeneracy at the Fermi 
level increases to 6-fold; in the presence of interactions 
the ground state can become triplet (5=1) for t! > 0.8. 
Hence we stop at t' < 0.8 where the ground state is 5 = 0. 

(ii) Interactions: As mentioned above, no transition 
occurs for U — Vij — 0. This further indicates that 
the transition is not a feature of the single-particle band- 
structure. 

(in) Commensur •ability: We have verified that the PEC 
state does not occur for electron densities different from 
p = |; for example, no transition to PEC or any "non- 
metallic" state occurs for 6 or 10 electrons on the 4x4 
lattice. 
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IV. APPLICATION TO REAL MATERIALS 

In the following we discuss how our theory applies to 
real materials, and may even give insight to the mech- 
anism of correlated-electron superconductivity in p = | 
materials. 



A. Application to organic CTS 

The superconducting organic CTS share many char- 
acteristics of other strongly-correlated superconductors, 
in particular the high-T c cuprates, including reduced di- 
mensionality and the presence of AFM near SC. At the 
same time, SC in the CTS occurs under pressure at a 
constant carrier density of p = | rather than under the 
influence of doping. A variety of exotic insulating states 
in addition to AFM, including CO 6 - and spin-gapped 
state o 66 ' 67 , as well as possible QSL states 68 , are proxi- 
mate to the superconducting state in the CTS. Our work 
shows that only the AFM phase is described by the dimer 
Mott-Hubbard model. Since in all cases the experimental 
systems are structurally related, with identical or near- 
identical molecular components, we believe that the same 
mechanism of SC should apply to them. In Reference 
we had pointed out how the PEC concept can perhaps 
lead to such a unifying theory. Here we expand on this 
theme. 



1. CTS with PEC insulating states 

We briefly review here experimental evidence for PEC 
formation in several 2D CTS families. 

(i) 6-(ET) 2 X: CO corresponding to the PEC and 
spin gap are found in the (9-(ET) 2 MM'(SCN) 4 family 66 . 
In M Af'=RbZn, the CO occurs below the MI transi- 
tion at T ~190K, while the SC appears below 20K 66 , 
The charge order pattern in the CO phase below the 
MI transition has been experimentally determined for 
MM'=RbZn and follows a horizontal stripe pattern (see 
Fig. 6 in Reference [60h . The horizontal CO is definitely 
not the WC. Rather, the CO pattern is precisely as ex- 
pected in the PEC, with • • • 1100- • • CO along the two 
directions of largest hopping (the p-directions in the 9- 
(ET)2X lattice), and ■ • • 1010- • • order along the direction 
of weakest hopping (c-direction in (9-(ET) 2 X). Experi- 
ments have revealed that with decreasing temperature, 
the c-axis lattice parameter decreases 7 ^. The decrease in 
the lattice parameter implies increased carrier hopping in 
this direction and therefore increased frustration within 
our theory, giving the transition to the singlet PEC and 
observed spin gap. 

(ii) a-(ET) 2 X: The crystal structure of a-(ET) 2 X is 
quite similar to that of #-(ET) 2 X. The existence of a SG 
opening below 136K has been known for some time in 
a-(ET) 2 I;A Below the 136K transition CO is found 



which has been confirmed to be of the same pattern as 
in 0-(ET) 2 X; see for example Fig. 2 in Reference l72l. 

(Hi) j3-(meso-DMET) 2 PF 6 : This 2D CTS exhibits a 
pressure-induced transition from CO to S C 73 ' 74 . While 
the charge order pattern in this CTS is referred to as 
"checkerboard" by the authors, the checkerboard pattern 
refers to meso-DMET dimers as units. In terms of meso- 
DMET monomer units, the CO pattern is the same as 
the PEC in Fig. H(c), with • • • 1100- • • in two directions 
and • • • 1010- • • along the third direction (see Fig. 2 in 
Reference [l!) 

(iv) f3' ' -X[Pd(dmit) 2 ] 2 : In this family the materials 
dmit molecules are arranged in dimers. The frustration 
varies with the cation X, with the least frustrated in the 
series showing AFM order—. Among the materials with 
larger frustration, X — EtMe^P has a SG transition at 
25K to what has been described as a VBS state 7 - 6 -. The 
experimentally determined bond and charge distortion 
patterns (see Fig. 3(b) in Reference [76|) are exactly as 
expected for the PEC, with period 4 charge and bond 
distortions along what is the x axis in Fig. [5] The in- 
tradimer charge disproportionation in particular argues 
against the SG state from being a simple VBS, which 
would require equal charge densities on all the molecules. 

2. K-(ET) 2 Cm(CN) 3 and EtMe 3 Sb[Pd(dmit) 2 ]2: QSL or 
charge- disproportionated states? 

There has been much recent interest in these CTS^ 8 - 
specifically because they present possible realizations of 
the long awaited QSl4£~— . In both cases the materi- 
als have nearly isotropic triangular lattices of dimer unit 
cells (corresponding to Fig. HJb) with t' = 1) within an 
effective p = 1 model. In k-(ET) 2 Cu 2 (CN) 3 (hereafter 
k-CN) the estimate for the Heisenberg exchange integral 
between n.n. dimers is J -220-250 K— i 7 - 8 -. ^ NMR 
experiments find absence of long range magnetic order 
down to 32 mK — . Very similar behavior is also seen in 
EtMe 3 Sb[Pd(dmit) 2 ] 2 (hereafter dmit-Sb) 7 ^. The insu- 
lating ground state in k-CN is close to SC, transition to 
superconductivity occurring under moderate pressure^ -. 
SC is found in the dmit family as well 8 - 1 -. 

Recent experiments in both k-CN and dmit-Sb have 
found peculiarities that appear to be unexpected within 
QSL theories. Below we list experiments that seem to 
indicate that apparent QSL behavior at low tempera- 
tures is giving way to a "hidden order", and perhaps 
even charge disproportionation, which by itself would be 
against spin-only models. 

(i) A second order phase transition is seen at 6 K in k- 
CN in measurements of heat-capacity Cjr^, 13 C NMR re- 
laxation rate 1/Ti— , and lattice expansion coefficients 8 - 4 -. 
The last experiment finds strong lattice effects at the 
transition, indicating possible role of charge degrees of 
freedom 84 . A symmetry-breaking and/or topological or- 
dering transition at T < 1 K has also been observed in 
dmit-SbS. 
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(ii) The specific heat C p in k-CN is linear in T for 
T between 0.75 and 2.5 K, indicating a gapless energy 
spectrum 8 - 2 -. The Sommerfeld coefficient 7 is nonzero and 
large even at T —75 mK 82 . Equally perplexingly, C p 
is independent of magnetic field up to 8 T, indicating 
absence of Zeeman coupling of spins to the field 8 ^. In 
contrast, thermal conductivity measurements down to 80 
mK indicate a spin gap^ 5 -. 

(iii) The temperature dependence of the thermal con- 
ductivity of dmit-Sb suggests gapless mobile excitations, 
but magnetic field dependence of the thermal conductiv- 
ity again indicates a gap£&. Taken together, (ii) and (iii) 
suggest gapless spin-singlet excitations but gapped spin- 
triplet excitations in re-CN and perhaps also dmit-Sb. It 
is conceivable that the gap in dmit-Sb is nodal 8 -^. 

(iv) Measurements of dielectric response in k-CN have 
shown increasing and frequency-dependent dielectric con- 
stant below 60 K, and possible antiferroelectric order- 
ing of dipoles at T c ~6 K 8 -^. The latter requires un- 
equal site charges on the molecules within the dimer unit 
cell a 87 i 88 . It may be relevant in this context that 13 C- 
NMR experiments on both k-CN 8 - 3 - and dmit-Sb^ 9 . find 
unusual line broadenings at low T that cannot be as- 
cribed to disorder—, which might also indicate charge 
disproportionation. Similar line broadening at low T in 
EtMe 3 P[(dmit) 2 ]2 occurs at the transition to the PECffi, 
which in turn gives way to superconductivity under 
pressure 8 -^. 

Although it is as yet not entirely clear whether or not 
existing spin-liquid models 4 ^— with appropriate mod- 
ifications can explain the above anomalies, it appears 
that for both re-CN and dmit-Sb, spin degrees of freedom 
alone cannot describe the low temperature (below 6K 
in k-CN and below IK in dmit-Sb) properties. Rather, 
any description of the ground state must involve charge 
as well as spin degrees of freedom. Assuming that the 
above experiments and their interpretations are correct, 
the proper theoretical model should have built-in signif- 
icant electron-lattice coupling and should lead to charge 
disproportionation and excitation energy spectrum with 
gapless singlet excitations and gapped spin excitations. 
We believe that the highly frustrated p = \ model sat- 
isfies all of the above criteria, in addition to providing 
the starting point for a theory of superconductivity in 
the CTS (see below). T- linear specific heat but gapped 
magnetic susceptibility were noted in nonmetallic vana- 
dium bronzes 9 - 1 - as far back as 1978, and at the time 
was considered to be a distinctive proof for bipolarons. 
The nonzero spin-singlet degeneracy within this model 
comes from tunneling motion of the bipolarons, causing 
them to "flip flop" between equivalent configurations^. 
We propose that a similar mechanism is at play in the 
present case at large t' . Particularly in k-CN the lattice 
structure and the orientations of molecules are such that 
the spin-singlet bonds, which occur between monomers 
belonging to two neighboring dimers can flip flop be- 
tween the monomers (see Fig. 4 in Reference [D.) Within 
our proposed picture for k-CN, the transition at 6K is 



to the spin-singlet charge-disproportionated state, with 
short range fluctuating order. Magnetic excitations even 
in this state, however, require breaking the spin-singlet 
bonds. Note that the observed difference between k-CN 
and dmit-Sb is expected within the p = \ model, since 
the simple description as a triangular lattice of dimers is 
no longer enough and the detailed couplings between the 
monomers in the materials are indeed different because 
of their different crystal structures. 



3. Consequence of stronger frustration — paired electron 
liquid and superconductivity 

As mentioned above, numerical results for t' > 0.8 are 
not useful due to the highly degenerate ground state be- 
coming spin-triplet in our finite clusters. We have sug- 
gested elsewhere that the occupied spin-singlet bonded 
'1-1' sites of the PEC can be thought of as effective sin- 
gle sites doubly occupied by charge carriers, and simi- 
larly the pairs of '0-0' vacancies can be thought of as 
single vacant sites 9 - 2 -. Such a mapping would transform 
the PEC to an effective checkerboard CO with alternate 
sites (in the square lattice representation) occupied by 
double occupancies and vacancies. The effective Hamil- 
tonian that describes the checkerboard CO in this case 
is a p = 1 extended Hubbard model with weak attractive 
U whose origin is the exchange interaction that stabilizes 
the singlet bond within the original p — | Hamiltonian. 
The n.n. interaction within the effective Hamiltonian re- 
mains repulsive to simulate the checkerboard ordering of 
the effective doubly occupied sites. 

While such a mapping is not rigorous, similar map- 
ping of n.n. spin-bonded sites to double occupancies 
has been routinely used in the literature on bipolaron 
model o 93 ' 94 . The difference between our work and tradi- 
tional bipolaron models is that the spin-singlet bonding 
within our work is driven primarily by AFM correlations, 
while within the bipolaron models it is a consequence of 
effective attraction due to the overscreening of the e-e 
repulsion by e-p interaction ;) 94 ! 95 . We have investigated 
the consequences of stronger frustration 9 -^ in the P = \ 
anisotropic triangular lattice within the effective p = 1 
extended Hubbard model with attractive U— . We found 
a CO-to-SC transition within the effective model 9 - 2 -, sug- 
gesting a PEC-to-SC transition within the p = h model 
with repulsive interactions. Although this result is not 
a proof of transition to SC within the repulsive p = | 
model, it nevertheless is instructive and provides the di- 
rection for future research. 

Demonstration of SC within the actual model will 
require further work. What is interesting though is 
that the proposed scenario can give a unified approach 
to SC in all organic CTS, irrespective of whether the 
insulating state proximate to SC is AFM2£, CO^^ 
or VBS 81 . Recall that within existing mean field 
theories the AFM-to-SC transition is driven by spin 
fluctuations^—, while the CO-to-SC transition is driven 
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by charge fluctuations^. Even if we ignore that recent 
precise numerical calculations^— have demonstrated 
the absence of SC within the proposed spin-fluctuation 
models in this context (thereby raising doubts also about 
mean-field theory of charge-fluctuation mediated super- 
conductivity), different mechanisms of SC for struc- 
turally related materials with identical or near-identical 
molecular components appear to be unrealistic. Within 
our proposed scenario, the charge-ordered or VBS sys- 
tems that exhibit SC are PECs in the semiconducting 
state, while the AFM systems under pressure transform 
to PECs first and then to SC (although due to actual 
crystal structures the "width" of the PEC region could 
be narrow to vanishing within the latter). The SC phase 
within this scenario is a paired-electron liquid. There is 
considerable overlap between these ideas and the one pro- 
posed by Moulopoulos and Ashcroft for the continuous 
electron ga s 45 ' 46 . 

B. Application to inorganic p = \ materials 

The thrust of our work has been to understand within 
a unified theoretical approach the variety of exotic insu- 
lating states that are proximate to the superconducting 
state in 2:1 cationic or 1:2 anionic 2D CTS. We have 
shown here that the peculiarities of these materials orig- 
inate from the unique behavior of p = i in the presence 
of both strong e-e interactions and lattice frustration. It 
is conceivable that the complex behavior of apparently 
unrelated inorganic families can be understood within 
the same broad theoretical approach. We list two such 
classes of materials below to point out the hitherto un- 
noticed similarities between them and p — \ CTS. 

1. Layered cobaltates 

Layered cobaltates NazCoO^ have attracted wide at- 
tention because of their 2D structure, tunable carrier con- 
centration, and the occurrence of SC 3 - 9 -. The Co ions form 
a 2D triangular lattice, are in their low-spin state, and 
their valence ranges from Co 3+ at x = 1 to Co 4+ at 
x = 0. Charge carriers are S — \ holes on the Co 4+ 
sites and the hole density p = 1 — x. Trigonal distortion 
splits the occupied ti g orbitals into e' g and a\ g orbitals. 
LDA calculations have suggested that although for large 
x the e' g orbitals occur below the Fermi level and an a\ g - 
only description is valid, this description breaks down at 
small x where e' g orbitals can be nearly degenerat e) 99 ' 100 . 
In contrast, correlated-electron calculations find that the 
dig - e' g energy separation is positive and relatively large 
for all a j 101 ' 102 , suggesting that low energy excitations can 
likely be described within ai s -only single-band models. 

The temperature dependent magnetic susceptibility 
x{T) shows a peculiar p-dependence within the family, 
with small p (large x) exhibiting strongly correlated be- 
havior and large p (small x) exhibiting weakly correlated 



behavior-iS 3 .. We have recently shown that p-dependent 
x(T), exactly as seen in the cobaltates, is expected within 
the single-band extended Hubbard model on a triangu- 
lar lattice 1 ^ 4 -. Equally interestingly, x(^) behavior in 
Na x Co02 is very similar to that in the family of CTS as a 
whole, where also x(^) shows a systematic p-dependence 
that is understood within the single-band extended Hub- 
bard modeliS^. It is tempting to compare now the super- 
conducting states in hydrated cobaltates and CTS with 
this apparent similarity in mind. 

Although superconductivity in Na x Co02 • 2/H2O oc- 
curs at x ~ 0.35^, it is now established that the Co-ion 
valency here is determined not only by the Na content, 
but also by H3O 4 " ions. There have been several reports 
that SC here occurs over a very narrow range of hole 
density p, and that maximum T c occurs at or very close 
to Co-ion valency 3.5+, corresponding to p = I22£r22£. 
If this is confirmed from future experimental work, it 
would appear that like the CTS, cobaltates are yet an- 
other example of a frustrated 2D p = | superconductor, 
suggesting that the mechanism of SC in the two families 
is related. 



2. p = \ spinels 

The B sublattice in spinel compounds AB2X4 form a 
frustrated three-dimensional (3D) pyrochlore lattice and 
usually consist of transition metal cations that possess 
partially filled ti g <i-orbitals. For integer occupancies of 
d-electrons per B cation the geometrical degeneracy of 
the underlying lattice is often lifted by orbital ordering 
(00), leading to formation of spin-singlet dimers. Only 
four of the many spinel compounds are superconduct- 
ing, of which three have effective carrier density p = ^: 
LiTi204, CuRh2S4, and CuRh2Se4. In LiTi204 there 
is one d-electron per two Ti 35+ ions; in CURI12S4 and 
CuRli2Se4 the Rh 35+ ions have average d-hole occupancy 
of \. Jahn- Teller instability or OO can lead to occu- 
pancy of the same t2 g orbitals, making the filled bands 
exactly ^-filled, as in the superconducting organic CTS 
(the apparently different spinel superconductor CUV2S4, 
also has noninteger number of ci-electrons per 

V 3.5+ 

ion.) 

Yet another similarity between the organics and the 
P = \ spinels is the proximity of the superconducting 
state to exotic semiconducting states. Thus CUR2S4 and 
LiRh 2 04, both isoelectronic with CuRh 2 S4, undergo MI 
transitions that are accompanied by CO. In CUR2S4 the 
Ir ions are charge-ordered as Ir 3 +-Ir 3+ -Ir 4+ -Ir 4+ along 
specific directions 1 ^ 9 -. This would correspond exactly to 
the ■ ■ • 0011- • • CO in our notation here. As in the 2D 
PEC, these inorganic 3D systems are spin- gapped due to 
the formation of Ir 4+ -Ir 4+ singlet bonds. It is interest- 
ing to recall that similar singlet bonds between Ti 3+ ions 
had been propose d 91 ' 110 many years back within bipo- 
laron theories of Ti407 and LiTi204. As explicitly shown 
in our work here, the singlet bond formation is a nat- 
ural consequence of frustration and ^-filling. Work is 
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currently in progress to extend the PEC concept to the 
checkerboard lattice, thought to be 2D equivalents of the 
pyrochlore lattice, multiple orbitals per site 1 *^. 

V. CONCLUSIONS 

In summary, there is an extraordinarily strong ten- 
dency to form spin-singlets in systems with charge carrier 
concentration precisely |. Naturally, in p = ^ this spin- 
singlet state is accompanied by CO. The stability of the 
PEC derives from the commensurability of the PEC at 
P = \- In the anisotropic triangular lattice, the PEC con- 
sists of charge arrangements • • • 1100- • • in two directions 



and • • • 1010- • • in the third direction. Thus although 
p = i in principle is incommensurate on the triangular 
lattice, "separately commensurate" periodic charge ar- 
rangements are nevertheless possible. 
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